Microscopic model of superconductivity in carbon nanotubes 
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We propose the model of a manifold of one-dimensional interacting electron systems to account 
for the superconductivity observed in ropes of nanotubes. We rely on the strong suppression of 
single-particle hopping between neighboring nanotubes in a disordered rope and conclude that the 
tunnelling takes place in pairs of electrons, which are formed within each nanotube due to the exis- 
tence of large superconducting correlations. Our estimate of the transition temperature is consistent 
with the values that have been measured experimentally in ropes with about 100 metallic nanotubes. 
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Carbon nanotubes have nowadays a great potential for 
technological applications in devices at the nanometer 
scale. This makes very important the precise knowledge 
of their electronic properties, taken as individual struc- 
tures as well as when packed in the form of ropes. The 
first step of discerning their metallic or insulating be- 
havior has been completed, as the results obtained from 
theoretical considerations have been confirmed exper- 
imentally It turns out that the low-energy spectrum 
of a carbon nanotube may be gapless or not depending 
on the particular wrapping of the graphene sheet to form 
the tubular structure. 

The electron interactions are also known to modify 
significantly the transport properties of the nanotubes 
|^,||. In the case of metallic single- walled nanotubes, 
the one-dimensional (ID) character of the system leads 
to a strong correlation among the electrons, making 
the so-called Luttinger liquid 1^-^ the most appropri- 
ate paradigm to describe the electronic properties [p|-p| . 
Evidence of Luttinger liquid behavior has been found in 
measurements of the tunnelling conductance in ropes |^] 
and in individual single- walled nanotubes Q. 

Experiments have been also carried out to probe super- 
conducting correlations in the carbon nanotubes |]l2| , |l3| . 
By suspending them between superconducting contacts, 
large supercurrents have been measured in samples al- 
most made of a single-walled nanotube and in ropes 
1^^ . More recently, superconducting properties have 
been measured in ropes suspended between nonsuper- 
conducting, good metallic contacts jlj] . Evidence of su- 
perconducting fluctuations has been also obtained in in- 
dividual single- walled nanotubes jist . In the experiments 
presented in Ref. jl^ , a drop by two orders of magnitude 
in the resistance has been found, down to the minimum 
value consistent with the number of conducting channels 
in the rope. This feature, together with its suppression 
under a suitably high magnetic field, shows the existence 
of superconductivity inherent to the ropes of nanotubes 
. The value of the transition temperature varies from 
one sample to other, being below 1 K in the two where 
the resistance drop has been measured. 



In this Letter we develop a model to account for the su- 
perconductivity intrinsic to the ropes of nanotubes. The 
problem that faces any model trying to describe such ef- 
fect is twofold. In the first place, it is known that the 
single-particle hopping between neighboring nanotubes 
in a rope is strongly suppressed ||l^. This is because, in 
general, the different helical structure of the nanotubes 
leads to the misalignement of their lattices and, therefore, 
to the difficulty of conserving momentum for an electron 
hopping from one nanotube to the other. On the other 
hand, the effect of superconductivity cannot rely exclu- 
sively on the properties of the individual nanotubes, since 
any correlation in a ID system can only develop a diver- 
gence at zero temperature Q. 

The experiments on the proximity effect of Ref. p^ ] 
show anyhow the existence of sensible superconducting 
fiuctuations in carbon nanotubes. The precise value of 
the critical supercurrents found there can be only ex- 
plained by the presence of a short-range attractive inter- 
action coming from the coupling to the elastic modes of 
the nanotubc |l^]. The point is that, as long as Cooper 
pairs arc formed locally with zero total momentum, they 
may overcome the special difficulty that single electrons 
find to tunnel between neighboring nanotubes. 

According to the work of Ref. |jl^, the amplitude It 
for tunnelling between nanotubes with the same chiral- 
ity and orientation is of the order tx ^ 0.01 eV. In 
a disordered rope, the misalignement of the lattices of 
neighboring nanotubes introduces in general an addi- 
tional suppression of the hopping amplitude by a factor 
^ exp(— i?ao((5fc)^/4), where R is the radius of the tube, 
5k is related to the mismatch of the Fermi points, and ao 
is a parameter of the order ~ 0.5 A [|6|. For a typical 
nanotube radius i? = 7 A, this factor is ^ 0.005. 

The tunnelling amplitudes have to be compared with 
the energy scale at which the metallic nanotubes behave 
as ID objects. This is the scale below which the 
gapless linear subbands dominate the physical properties, 
and it can be estimated as £'c ~ 0.1 eV. The probability 
A2 of tunnelling of a Cooper pair, given by the square 
of It in units of i?c, has a relative weight of the order 
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~ 0.01. On the other hand, the probabihty of tunnelhng 
of a single electron has a relative weight of the order 
^ 0.0005. This means that pair hopping is the dominant 
process for tunnelling between neighboring nanotubes in 
a disordered rope. 

Prior to considering pair hopping, the rope can be de- 
scribed at low energies by a model in which each nan- 
otube is treated as a ID system, whose charge can inter- 
act with the charge in the other nanotubes of the rope. 
The hamiltonian for this model, including a collection of 
metallic nanotubes a = 1 , . . . n with linear branches of 
different chirality i = +,— , can be written in terms of 
the respective density operators paia [fSl : 
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Hi = -Vp / dfc^ : Paia(k)pata{~k) I 
aia 

+ 2/ dk^Pa^a{k)^Vab{k) pbja'{-k) (1) 
'^^ aia bjcr' 

where kc — Ec/vp. 

For the potential between different nanotubes Vab we 
take the Coulomb interaction Vc{k) — e^/(47r^) log |(fcc-l- 
k)/k\ p9| , which remains long-ranged in one spatial di- 
mension [|2^. According to the results of Ref. [ p7[ , 
the interaction potential Vaa within each nanotube in- 
cludes moreover the effective short-range attraction com- 
ing from the coupling to the elastic modes, so that 

Ka(fc) = eV(4^') log |(fcc + fc)/fc| -5/(2^)- The strength 
g of the attractive interaction is inversely proportional 
to the mass of the atoms and directly proportional to 
{t' /vsY, where t' is the modulation of the lattice hop- 
ping and Vs is the speed of sound . A rough estimate 
for a carbon nanotube gives g/vp ~ 0.2 . 

Terms which couple the spin densities have been 
neglected in writing the hamiltonian (|]). These are 
backscattering (BS) interactions, which arise as a rem- 
nant of the Coulomb interaction at short distances. Be- 
ing local in the nanotube lattice, their couplings are re- 
duced by a relative factor of the order ^ O.la/i?, where 
a is the nanotube lattice spacing ||9|-[ll|| . These terms are 
marginally relevant in the renormalization group sense. 
This means that they have greater strength as the model 
is scaled to smaller energies, but the rate of increase 
starts being quadratic in the own couplings. Thus, the 
theory has to be scaled to extremely low energies, many 
orders of magnitude below Ec, to have the BS couplings 
comparable to the couplings in JlOt . 

The bundle of ID electron systems coupled only by 
charge interactions ressembles the system proposed in 
Ref. for the description of the sliding Luttinger liq- 
uid. Here the Coulomb interaction couples each of the 
nanotubes to all the others in the bundle. We recall 
that, for the samples of the experiments reported in Ref. 
p^ , the number of metallic nanotubes is very large, of 
the order n ~ 100. Then, the coupling of the charge in 



the different nanotubes leads to a significant reduction 
of the effective repulsive interaction. On the other hand, 
the effect of the intra-tube attractive interaction does 
not depend on the number of nanotubes and, in certain 
regime, it may dominate over the Coulomb interaction. 

The above point can be checked by looking at the cor- 
relators of the model governed by Hi. These can be 
computed exactly by means of bosonization techniques. 
For instance, the propagator D^^^ {x, t) for the Cooper 
pairs within each nanotube factorizes into the different 
interaction channels, taking the form 

Z?(")(x,t) ^ (vI/+^^(x,t)vI/+_^(:r,i)*,+T(0,0)*,_^(0,0)) 
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C{x,t)\{N{x,t)\{F{x,t) 



(2) 



The first factor in Eq. (g) corresponds to the contri- 
bution of the total charge density, which is given at zero 
temperature by the expression 

C{x, t) = exp I I dk — (1 — cos(fc2:) cos{vpkt)) 

y 2n Jq p[k) k 

(3) 



where fi{k) — 1/ \J\ + %nVc(k')l'^F — '^g/i'T^vp) and 
vp{k) = vp/p{k). The factors N{x,t) correspond to the 
rest of the charge channels, and they have a form simi- 
lar to (^) but with p = — ■ig/{'KVp) instead of p{k) 
andvp ~ vp/p instead of vp{k). The factors -F (a;, i) cor- 
respond to the noninteracting channels with p — 1. The 
computation can be extended to the case of temperature 
T 7^ 0, just by inserting the factor l+2/[exp(u_F|fc|/T)-l] 
in the integrand of expressions such as (|^) ||] . 

The couplings that produce the most accurate fit 
of the critical supercurrents reported in Ref. |^2| are 
2e^/{Tr^vp) « 1.0 and Ag/invp) w 0.6 0. For n - 100, 
the strong reduction of the Coulomb interaction implies 
that the effect of the intra-tube attractive interaction pre- 
vails in the system. In the g-ology description Q, the 
model is in the regime with short-range attractive cou- 
pling g2 < 0, where the singlet superconductor response 
is enhanced over the charge-density-wave response. The 
temperature dependence of the Fourier transform D^*^^ of 
the propagator at zero frequency and momentum is rep- 
resented in Fig. 0. The enhancement of the propagator 
is the signal that large superconducting correlations exist 
in the individual nanotubes at low temperatures. 

The above considerations are pertinent to the system 
at half-filling, in which the linear subbands cross at the 
Fermi level. When the nanotubes are slightly doped, the 
shift of the Fermi level gives rise to four Fermi points. 
The Cooper pairs have then the possibility to resonate 
between the outer and the inner gapless subbands. Any- 
how, as long as the repulsive interactions mediating these 
processes are reduced by a relative factor as small as that 
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of the BS interactions |^,|0[ , the s-wave pairing is favored 
over other channels with different symmetry. 
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propagator of the Cooper pair from a metaUic nanotube 
a to another 6 as a function of the distance x along the 
rope and the positions la and 4 of the nanotubes in the 
transverse section of the rope, D{la, W, x, t). This object 
can be related to the propagator I^^^^ along each nan- 
otube through the self-consistent diagrammatic equation 
in Fig. |2| which takes into account the dominant terms 
in powers of A2 . 
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FIG. 1. Plot of the propagator J)*-"' at zero frequency and 
momentum versus T / Ec, for 2e^ / {■k'^vf) ~ 1.0 . The dashed 
line corresponds to the case n = 1 and g = Q, and the solid 
lines to n = 100 and the respective values (from top to bot- 
tom) Ag/ijiVF) =0.75,0.5,0 . 

The tunnelling of Cooper pairs between the nanotubes 
can be taken into account by modifying the hamiltonian 
ffu) with the additional term 



b I dk dp dp' 



{a,b) 



ka 



[^+^ik + p)^+^^^{~p)^i,nik + p')^b-ni-p') + i^-^) (4) 



where ^^j^. is the electron operator and the sum runs over 
all pairs (a, 6) of nearest-neighbor metallic nanotubes. 

The term (j^ is a relevant perturbation from the renor- 
malization group point of view. However, the anomalous 
scaling dimensions of the relevant perturbations turn out 
to be in general rather small. They can be computed in 
the boson representation and, in the case of the term (||), 
the result is 72 = 2[l/{4n^i{ko)) + {n - l)/(4n/x) - 1/4]. 
ko represents some effective value, which does not affect 
significantly the estimate. For the above mentioned cou- 
plings, we obtain 72 « —0.07. Thus, even four orders of 
magnitude below Ec, we observe that A2 is not enhanced 
by more than a factor (lO^'*)''^ ^ 2. 

The proposed model has a superconducting instabil- 
ity at some finite temperature, provided that the Cooper 
pairs are able to percolate in the transverse directions of 
the rope. This is the case of the superconducting ropes of 
Ref. ]T^ , although to meet such experimental condition 
one has to find the appropriate sample out of a large num- 
ber of them [ p3[ . Assuming the coherence in the hopping 
of pairs in the transverse directions, we may write the 
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FIG. 2. Self-consistent diagrammatic equation for the 
propagator D of Cooper pairs along the rope. 

By introducing the Fourier transform with respect to 
{x,t) as well as in the la variables ||2^, the equation in 
Fig. ^ can be written for the Fourier transformed prop- 
agator D in the form 

D{q- k, ujk) = 5(0) (fc, uu) -f (fc, ujk)\2 {q)D{q- k, Uk) 

(5) 

We are interested in the propagation of the Cooper 
pairs from a metallic nanotube to the rest in the rope, 
which is given by D{0; 0, 0) = J dx J dt D{la,lb', x, t). 
By solving Eq. (0), we get the result 
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1-A2(0)i?(o)(0,0) 



(6) 



The relevant dependence of 13(0; 0,0), as well as of 
-D*^"^ (0, 0), is on the temperature. It becomes clear that, if 
the superconducting correlations are such that i:>(o)(0,0) 
has a divergence at T = 0, then the effect of pair hopping 
gives rise to the appearance of a pole at a finite value of 
T in the propagator of the Cooper pairs in the rope. 
According to the conventional interpretation, this is the 
signal of the condensation of Cooper pairs and the onset 
of the superconducting transition in the system. 

In practice, when dealing with a rope of finite length 
L, the divergence of l)(°)(0, 0) is cut off at a temperature 
scale about one order of magnitude below vp/L. The 
curves shown in Figs. ^ and ^, for instance, have been 
obtained for a system with L — 1000//cc, which corre- 
sponds to L ^ 1 /im with an appropriate choice of the 
length scale. 

Taking our estimate of the pair-hopping amplitude 
A2 ^ 0.01, we observe that the length L = lOOO/fcc may 
be in some samples at the limit below which a super- 
conducting instability cannot arise in the system. This 
depends on the spatial distribution of metallic nanotubes, 
which should determine more precisely the effective value 
of A2 to be used. Quite remarkably, a superconducting 
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transition has been found in samples whose length is 1 
/im or greater, while a sample with low resistance and 
0.3 jjm long has shown no transition at all jl^. 



should be studied since, as observed from Figs, ^and ||, a 
slight change in that parameter may result in an increase 
of Tc by more than one order of magnitude. 
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100 



D 



(0) 



10' 




10 



10 



10 



10 



FIG. 3. Logarithmic plot of at zero frequency and 

momentum versus T /E^- The solid lines are the transposition 
of the solid curves in Fig. ^ for n = 100 and the respective 
values Ag/{-KVF) = 0.75,0.5,0 . The dashed line corresponds 
to the case n = 400 and Ag/ljivp) — 0.75 . 

Relying also on our estimate for A2, we observe from 
Fig. H that the temperature of the transition to the su- 
perconducting phase in a disordered rope with about 100 
metallic nanotubes can be in the range between 
and 10~^£'c. As the natural energy scale in our model is 
Ec ~ 0.1 eV, this sets the scale for typical transition tem- 
peratures in the range between 0.1 K and 1 K . These 
values are consistent with the transition temperatures 
Tc « 0.1 K and Tc ~ 0.4 K measured experimentally in 
the two samples of Ref. Q . 

The results of our analysis show that the superconduc- 
tivity in the ropes of nanotubes is close in nature to that 
of the alkali-doped fuUerenes . We have seen that the 
tunnelling of electrons between neighboring nanotubes in 
a rope takes place in pairs, which are formed within each 
nanotube due to the large superconducting correlations 
which develop at low temperatures. 

The low values of Tc compared to those of the alkali- 
doped fuUerenes can be understood on phcnomcnological 
grounds by the fact that the elcctron-phonon coupling is 
smaller in the nanotubes. This is consistent with a higher 
estimate of Tc obtained in Ref. ||l^ from measurements 
on small-diameter single-walled nanotubes, which have 
a larger electron-phonon coupling than the nanotubes in 
the ropes. As seen in Fig. ||, another way to obtain a 
higher Tc in the ropes may be to increase the number 
of nanotubes of the samples. Finally, the feasibility of 
controlling the rate of tunnelling between the nanotubes 
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